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Abstract 

It has previously been proved that the lowest order supersymmetric WKB approx- 
imation reproduces the exact bound state spectrum of shape invariant potentials. We 
show that this is not true for a new, recently discovered class of shape invariant po- 
tentials and analyse the reasons underlying this breakdown of the usual proof. 
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In the past few years, the supersymmetry inspired WKB approximation (SWKB) 
[1] has received considerable attention [2]. One reason for this was the proof [3] 
that the leading order SWKB quantization condition reproduces the exact bound 
state spectra for any shape invariant potential (SIP) [4]. Subsequently, Adhikari et 
al.[5] showed by explicit calculation that the higher order corrections to the energy 
eigenvalue spectrum vanish to 0{h^) for all then-known shape invariant potentials 
[6]. For these cases, all the higher-order corrections have since been shown to vanish 
[7,8]. The SWKB quantization condition has also been applied to many non-shape 
invariant potentials and it turns out that even though the SWKB formula does better 
than the usual WKB approach in most cases [5,9,10], it has never been found to be 
exact for these, thus suggesting that perhaps shape invariance is not only sufficient but 
even necessary for the lowest order SWKB formula to yield exact energy eigenvalues 
[10]. Barclay and Maxwell [8] proposed a very simple condition on the superpotential 
which ensures that the lowest order SWKB formula gives exact bound state spectra 
with all higher order corrections zero. By analysing this simple condition, they have 
suggested that there are no SIPs other than those tabulated by Dutt et al.[6] and 
the one example found by Levai [11]. It should be noted here that most of these 
potentials are also contained in the hst of Infeld and Hull [12]. 

Recently, a large class of new SIPs has been discovered by us [13,14] disproving the 
conjecture that no more examples exist. These new SIPs are reflectionless and possess 
an infinite number of bound states. They can be looked upon as q-deformations of the 
symmetric Rosen-Morse potentials corresponding to one or multi- soliton solutions. 
Although these new potentials cannot be written in a closed form using elementary 
functions, the energy eigenvalue spectrum is known analytically. It is then reasonable 
to enquire if the lowest order SWKB formula yields the exact bound state spectra for 
these new SIPs or not. It is worth emphasizing that the answer to this question is 
not obvious. On the one hand, since the lowest order SWKB formula has been shown 
to be exact for all SIPs [3], it is natural to expect that even for the new SIPs, the 
lowest order SWKB formula must be exact. On the other hand, Barclay and Maxwell 
[8] have claimed that there are no other potentials for which the lowest order SWKB 
formula is exact and the higher order corrections are zero. Thus, it is clearly of 
interest to know if the lowest order SWKB formula is exact for the new SIPs or not. 
This is the question that we address in this letter. 

In SUSY quantum mechanics, the partner potentials V±{x, Qq) = ± W'{x, ao) 
are said to be shape invariant if (2m = 1) 

W^{x, ao) + nW'{x, ao) = W\x, ai) - hW'{x, ai) + R{ao), (1) 
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where W is the superpotential, oq denotes a set of parameters, oi = /(oo) is an 
arbitrary function of and -R(ao) is independent of x. For all the standard SIPs 
[6,11] - for which the lowest order SWKB formula is exact - it turns out that the 
parameters Oi and Oq arc related by translation (ai — ao + ah; a being a constant). 
On the other hand, for the new SIPs, ai and ao are related by scaling (ai = goo; 
< g < 1). Expansion of the superpotential and R{ao) in powers of ao gives 

oo oo 

W{x,ao) = J29j{x)al, R{ao) = ^Rja^, ai = qa^. (2) 

Substitution into the shape invariance condition, eq. (1), and equating powers of ag 
permits a full determination of the quantities gj{x) for any choice of the constants 
Rj. A particularly simple choice is to take only one non-zero Rj. It has been shown 
[13,14] that this special choice gives the self-similar refiectionless potentials studied 
by Shabat [15] and Spiridonov [16]. In that case [13,14] 

Ri 

Rn — Rl^nl: Pi — J^—, (3) 

(1 + g) 

= /3„x-- = ( - p„i\-;(f|^„) g'aA-,)x"-'. (4) 

so that W{x) is an odd function of x and hence SUSY is unbroken. The exact bound 
state spectrum for V-{x, ao) is now easily calculated and is given by [13,14] 

K^i^iaoV^' ^ = 0,1,... (5) 

For the case of unbroken SUSY, the lowest order quantization condition is 

J VE -W^ dx = mih, n = 0, 1,... (6) 

By construction this will always give the exact ground state energy {Eq — 0), but 
what about the excited states ? Since W{x) is not available in closed form in terms 
of elementary functions, the SWKB integral has to be done numerically. A method 
for calculating W{x) has been described before [14] and careful numerical integration 
shows that for all q e (0, 1) the lowest order SWKB formula is not exact for the new 
potentials. The left hand side of equation (6) can be thought of as defining a function 
n{E) = (nh)"^ J \/E — W^dx which will take on the integer value n ioi E — En, if 
this lowest order quantization condition is exact. In Table 1 we display the results 
of doing the integral using the exact results for E = Ei, i.e. calculating n{Ei). For 
q close to one the discrepancy is small (though still significantly greater than the 
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numerical uncertainties), but it diverges (roughly as In 5) as g — > 0. At exactly q — 0, 
the Rosen-Morse potential, the condition is known to be exact, a discontinuity that 
is presumably a reflection of the way in which the infinitely many excited states for 
q > collapse to a single one when q = 0. 

The limit g — > 1 is much smoother and corresponds to the harmonic oscillator, for 
which the SWKB formula is also known to be exact. In fact, it appears that one can 
perturb about this limit for small values of e = 1 — g. The shape invariance condition 
(1) can be expanded to find W as a power series in e and this used to evaluate the 
integral. One finds that 

r>'iEi) = l + ^e' + 0{e% (7) 

which accords well with the numerical results close toq — 1. All the above calculations 
are done for the simplest case where only Ri is taken non-zero. It is easily checked that 
the lowest order SWKB formula does not give correct eigenvalues for the more general 
and complicated cases of new SIPs corresponding to several non-zero Rj either. This 
is of course intuitively expected. 

What is the special feature of these new potentials that interferes with the proof 
that condition (6) is exact for shape invariant potentials ? To understand this, let us 
take a fresh look at the derivation of the lowest order SWKB quantization condition, 
paying particular attention to the role of h. The lowest order WKB quantization 
condition for the potential V_{x, Gq) is (2m = 1) 

- V_ {x, oo) dx^{n + l /2)7rh. (8) 

In terms of the superpotential W, this reads 

J - W^{x, ao) + hW'{x, ao) dx ^ {n + l/2)7ih. (9) 

Now comes the crucial step in the derivation. One argues that formally W'^ is of 
0{h^) while hW is of 0{h) and hence one can expand the integrand on the left 
hand side in powers of h to get the lowest order SWKB quantization condition (6). 
This assumption is justified for all the standard SIPs [6,11] since is indeed of 
0{h^) while hW' is indeed of 0{h). One might object to this procedure since the 
resulting potential V- is then /?.-dcpcndcnt. However, in every one of those cases this 
^-dependence can be absorbed into some dimensionful parameters in the problem. 
For example, consider 

W^Atanhx (10) 

so that 

V4x) = - A{A + h)sech^x. (11) 
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Taking A such that A{A + h) is independent of h gives the desired /i-independent 
potential (the additive constant is irrelevant and so can contain h). Such a move may 
appear to be of limited value since one cannot apply SWKB directly to a superpo- 
tential W which is now /i-depcndcnt. However, because A is a free parameter, one 
can continue the SWKB results obtained for A (and hence W) independent of h over 
to this superpotential and so obtain an SWKB approximation for a ^-independent 
potential. 

What about the new potentials ? In the simplest of these cases, the only free pa- 
rameter in the problem (apart from q) is the combination Riao, on which W depends 
as W{x, Ritto) = y/RiaoF{y/Riaox/h). Incorporating different dependences on h in 
RittQ will give different ones in W, and En, but F is a sufficiently complicated 
function that there is no way of eliminating h from W^. This is a direct consequence 
of the scaling reparametcrisation ai = qao not involving h: if W'^{x, oq) were indepen- 
dent of h, so would W'^{x, tti) be and in taking the lowest order of (1) one would get 
W^{x, ao) = W^{x, Oi), which corresponds to the harmonic oscillator. A reparametcr- 
isation of the form ai = ao + ah clearly avoids this. Furthermore, with ai = qao, 
and hW are now of a similar order in h. The basic distinction between them involved 
in deriving equation (6) is thus no longer valid and we are prevented from deriving 
the SWKB condition for these new potentials. Similar difficulties will arise for po- 
tentials with more Rj non-zero and for the other ^-independent reparameterisations 
considered in [14]. 

It might be thought that an alternative derivation of the SWKB series for these 
potentials could be constructed by considering W — y/R^F{y/R^x/ H), where H 
is simply a free constant unrelated to h. The series can be safely derived for all H since 
W is ^-independent and the result continued to H — h at the end. Unfortunately, 
when this is done it becomes clear that the structure of all the terms is such that 
the entire series is then all of 0{h), so that one cannot usefully separate out a lowest 
order condition corresponding to (6) from a set of higher order corrections. It is then 
little surprise that an arbitrarily selected piece (6) of the full condition does not give 
the exact spectrum. 

Since it appears to have been overlooked in the extensive literature (see e.g. [17]) 
on symmetric, refiectionless potentials containing infinitely many bound states, we 
point out that the inextricable /i-dependence of these new shape invariant potentials is 
actually a general consequence of them possessing these properties. Reflect ionlessness 
implies that a potential is finitely deep, and the infinite spectrum means that the top 
of it is given by the limit of as n — oo. However, as /i — * 0, so the 

potential must also distort in the semi-classical limit; a little further reflection shows 
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that it always tends towards a free-particle potential. A similar conclusion can be 
reached for the A^-soliton rcflectionless potential (which contains bound states) by 
explicitly restoring h ^ I'm the inverse scattering formalism for these potentials [17]. 

We have thus seen that the SWKB quantization condition (6) is not the correct 
lowest order formula in the case of the new SIPs and hence it is not really surprising 
that eq. (6) does not give the exact eigenvalues for these potentials. In other words, 
it remains true that the lowest order SWKB quantization condition is exact for SIPs 
(if the SUSY is unbroken), but only in those cases for which the formula is applicable 
in the first place. It is thus still the case that the SIPs given in refs.[6,ll] are the only 
known ones for which the lowest order SWKB formula is exact and the higher order 
corrections arc all zero. 

It is a pleasure to thank R. Dutt, A. Gangopadhyaya, C.J. Maxwell and A. Pag- 
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partment of Energy under grant DE-FGO2-84ER40173. 
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0.01 


1.238 


0.1 


1.101 


0.2 


1.063 


0.3 


1.042 


0.4 


1.029 


0.5 


1.019 


0.6 


1.012 


0.7 


1.007 


0.8 


1.003 


0.9 


1.001 



Table 1: Results of evaluating the lowest order SWKB integral using the exact value 
El — Ritto for different values of q. An exact quantisation condition would give n — 1 
in every case. Note that the RiUq dependence cancels in the integral, so these results 
do not depend on this free parameter. 
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